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Abstract 
O 

We analytically work out the long-term orbital perturbations induced by the leading 

order of perturbing potential arising from the local modification of the Newton's inverse 
bJ[)' square law due to a topology M. 2 x S 1 with a compactified dimension of radius R recently 

proposed by Floratos and Leontaris. We neither restrict to any specific spatial direction k 

p^i ' for the asymmetry axis nor to particular orbital configurations of the test particle. Thus, 

^ , our results are quite general. Non-vanishing long-term variations occur for all the usual 

l/"*) ' osculating Keplerian orbital elements, apart from the semimajor axis which is left unaffected. 

By using recent improvements in the determination of the orbital motion of Saturn from 
Cassini data, we preliminarily inferred R > 4 — 6 kau. As a complementary approach, 
the putative topological effects should be explicitly modeled and solved-for with a modified 
version of the ephemcrides dynamical models with which the same data sets should be 
reprocessed. 
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l j_! '. 1 Introduction 



In [T] non-trivial topological modifications of the Newton's law were proposed on the ground 
of certain considerations of astrophysical and cosmological nature pertaining various forms of 
anisotropics occurring at such scales. In particular, Floratos and Leontaris [T] looked at local 
changes of the topology of the Euclidean space from M 3 to M 2 x S 1 , with a compactified dimension 
with scale R. 

In Section [21 we will analytically work out the orbital effects induced by them on the motion 
of a test particle moving around a central body of mass M. For distances r smaller than the 
compactification radius R, the correction to the Newtonian potential £/n is [I] 
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1 Floratos and Leontaris 1 discussed also the case r S> R- 



where G is the Newtonian constant of gravitation, Q is the Riemann Zeta function, P^j is 
the Legendre polynomial of degree 2j, and cos 9 is the cosine of the angle between the unit 
position vector r of the point in space at which the potential is evaluated and the direction k 
of the asymmetry axis. Contrary to Floratos and Leontaris [I], we will neither align k with the 
reference z axis nor restrict to any specific orbital configuration for the test particle. 

Since it is expected that R is quite large, in the following we will consider only the effects 
induced by the leading order term (j = 1) in eq. ([1]). Such an assumption will be a-posteriori 
justified in Section [31 where we preliminarily compare our results to the latest observational 
determinations of the orbital motion of Saturn. 

Section U] is devoted to the conclusions. 

2 The long-term rates of change of the osculating Keplerian 
orbital elements 



The long-term rates of change of the usual osculating Keplerian orbital elements due to the 
leading order term of C/ per t (j = 1) can be straightforwardly worked out with standard pertur- 
bative techniques. For example, by averaging the leading order term of eq. ([1]) over one orbital 
revolution, the Lagrange equations for the variation of the elements [2] yield 
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where [2] a,e,I,D,uj,w = J7 + w, 7W are the semimajor axis, the eccentricity, the inclination to 
the reference {x, y} plane of the coordinate system adopted, the longitude of the ascending node, 
the argument of pericenter, the longitude of pericenter and the mean anomaly, respectively; 
nb = \JGMjc? is the Keplerian mean motion. The simplifying condition e, / ~ was not 
imposed in obtaining eq. ([2])-eq. (jlip . which are, thus, exact. 

According to eq. (|3|), for circular orbits (e = 0) the eccentricity is left unaffected. 

3 Confrontation with the observations 



By using eq. ([2])-eq. (jlip it is possible to infer lower bounds on R from solar system plane- 
tary orbital motions. Recently, Fienga et al. [3] determined supplementary precessions of the 
perihelia and the nodes of some planets of our solar system which, in principle, account for 
any unmodeled/mismodeled dynamical effects with respect to those accounted for in the force 
models fitted to the observations in producing the INPOPlOa ephemerides which use the mean 
terrestrial equator at J2000.0 as reference {x, y} plane. Since the predicted rates of change are 
proportional to a 3 ' 2 , the largest effects are sensed by the outer planets. Thus, we will use the 
perihelion and the node of Saturn, whose precessions are nowadays known at a ~ 0.4 — 0.6 mas 
cty^ 1 level of accuracy [3]. 

In Figure Q] we display the lower bounds on R, in kau, obtained from the Kronian perihelion 
and node as functions of the right ascension a and declination 5 of the asymmetry axis k. From 
an inspection of Figure [TJ it turns out that the perihelion of Saturn yields up to R > 6 kau for 
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Figure 1: Minimum compactification scale R from the supplementary advances of the perihelion 
and the node of Saturn by Fienga et al. [3] as a function of the right ascension a and the 
declination 5 of the asymmetry axis k. We adopted eq. ([9]) and eq. ([7]). 



most of locations in the sky of the asymmetry axis. The perihelion-inferred bounds on R are 
weaker in a few positions of k, but the node can be fruitfully used in a complementary way to 
infer tight constraints also there. Thus, we can reasonably conclude that R > 4 — 6 kau for all 
over the sky. 

4 Conclusions 

We looked at some potentially observable consequences of a possible non-trivial local modifi- 
cation M? x S 1 of the topology of the Euclidean space with a compactified spatial dimension 
with scale R recently proposed by Floratos and Leontaris. We analytically worked out the or- 
bital effects due to the first term in the expansion of the resulting perturbing potential. Apart 
from the semimajor axis, all the other standard osculating Keplerian orbital elements undergo 
non-vanishing long-term variations. We did not make a-priori simplifying assumptions on the 
orbital configuration of the perturbed test particle. Moreover, we did not choose any specific 
spatial direction for the anisotropy axis. 

We preliminarily inferred lower bounds on R from the latest improvements in the determi- 
nation of the orbital motion of Saturn from Cassini data. We obtained R > 4—6 kau, depending 
on the orientation of the asymmetry axis in the sky. With such figures it is possible, a pos- 
teriori, to justify our initial choice of considering just the leading order term in the expansion 
of the perturbing potential since the precessions of higher order, computed for them, would be 
far smaller that the present-day level of accuracy in determining the extra-rates of the node 
and the perihelia of Saturn. We did not consider other known mismodeled/unmodeled forces 
causing potentially competing precessions like the Sun's oblateness, the PPN /3 parameter and 
the Lense-Thirring effect since, at present, they are negligible for Saturn with respect to the 
current level of accuracy in constraining its supplementary precessions. 

For a more refined analysis, it would be possible, in principle, to employ an ad-hoc modified 
version of the usual suite of standard force models to be fitted to the same data records, and 
explicitly solve for R or some other equivalent adjustable free parameter. 

As a possible further investigation of the viability of the suggested modified topology, the 
behavior of a typical Oort comet may be studied in order to see if the usual picture of the 
Oort cloud, which extends well beyond the previously inferred lower bound for R, would be 
notably altered. 
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